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The head-on collision of identical neutron stars from rest at infinity requires a numerical simulation 
in full general relativity for a complete solution. Undaunted, we provide a relativistic, analytic 
argument to suggest that during the collision, sufficient thermal pressure is always generated to 
support the hot remnant in quasi-static stable equilibrium against collapse prior to slow cooling via 
neutrino emission. Our conclusion is independent of the total mass of the progenitors and holds 
even if the remnant greatly exceeds the maximum mass of a cold neutron star. 
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I. INTRODUCTION 

Binary neutron stars are among the most promis- 
ing sources for gravitational wave laser interferometers 
now under construction, like LIGO, VIRGO, GEO and 
TAMA. This fact has motivated an intense theoretical ef- 
fort to understand their dynamical evolution and predict 
the gravitational waveforms emitted during their inspiral 
and coalescence. 

Here we focus on the final collision and coalescence of 
the two stars. The complex interplay between hydrody- 
namics and gravitation, even in the idealized adiabatic 
case, requires that this late epoch be tackled by numer- 
ical means. The numerical integrations must be carried 
out in full general relativity even for a qualitative, let 
alone a quantitative, understanding of the key phases of 
the final evolution. Strong-field phenomena characteriz- 
ing the final stages of binary neutron star evolution, like 
gravitational collapse leading to black hole formation, are 
not addressed at all by Newtonian theory. 

Consider the following three strong-field questions: Do 
binary neutron stars collapse to black holes prior to con- 
tact? If not, does the remnant undergo gravitational 
collapse immediately following the merger whenever its 
rest mass exceeds the maximum equilibrium mass of an 
isolated cold, rotating neutron star, or is the collapse de- 
layed ? What is the final fate of a collapsing remnant if its 
angular momentum exceeds the Kerr limit, J/M 2 < 1 ? 
None of these questions can be answered via a Newtonian 
simulation. The answers have important consequences 
for the design of advanced detectors and theoretical tem- 
plates for measuring gravitational waves during the late 
epochs following collision, as well as other astrophysical 
consequences (e.g., gamma-ray burst models). 

Here we address the second of the three questions 
posed above. To simplify the analysis we treat the sim- 
plest merger scenario - the head-on collision of two iden- 
tical, nonspinning, cold neutron stars which fall together 
from rest at infinity. The stars may have a combined rest 
mass which greatly exceeds the hydrostatic equilibrium 
mass limit for an isolated, nonrotating, spherical star 
governed by the same cold equation of state. Even this 
problem requires a numerical simulation in full general 
relativity for a complete solution. However, by making a 
suitable set of simplifying approximations, we present an 



analytic argument to determine the fate of the remnant 
immediately following merger. Our calculation may pro- 
vide a useful comparison for a 2D or 3D relativistic code 
which can perform such a simulation. 



II. ALTERNATIVE SCENARIOS 

Following the head-on collision of the two neutron 
stars, two recoil shocks form at the point of contact and 
propagate back into each star along the collision axis. 
Matter passing through the shock is heated and decel- 
erated, as its kinetic energy is converted into thermal 
energy. The configuration is not in equilibrium following 
the passage of the initial recoil shocks, but undergoes os- 
cillations, damped by acoustic (viscous) dissipation and 
gravitational radiation reaction forces. Only a small frac- 
tion of the total rest mass escapes from the system and an 
even smaller fraction of the total mass-energy is carried 
off by gravitational waves. 

Two questions arise: Is the thermal pressure generated 
by shock and acoustic dissipation sufficient to support 
the merged remnant in hydrostatic equilibrium prior to 
cooling by neutrino dissipation? Does the answer depend 
on the total rest mass and its relation to the maximum 
rest mass of an isolated, equilibrium mass supported by 
the same initially cold equation of state? 

Consider two alternative scenarios. In the first we sup- 
pose that the shock-generated thermal pressure is suf- 
ficiently large to prevent the remnant from undergoing 
immediate collapse, allowing it instead to settle down to 
quasi-equilibrium in a few dynamical timescales (~ msec) 
following the initial oscillations. By then, the bulk of the 
kinetic energy will have been converted into thermal en- 
ergy, the rest having been emitted as gravitational waves. 
The star will relax to a hot, quasi-static, spherical equi- 
librium state. This equilibrium state will last many dy- 
namical timescales because the thermal neutrinos, which 
eventually carry off the thermal energy, leak out slowly, 
on a neutrino diffusion timescale (~ 10 sec). As the star 
cools, it slowly contracts. If its total rest mass is less than 
the critical maximum value, it will evolve to a static, cold 
endpoint configuration and cease all further contraction. 
If its rest mass exceeds the critical value it must ulti- 
mately undergo catastrophic collapse to a Schwarzschild 
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black hole. Since the quasi-static configuration is spheri- 
cal and not too far from the Schwarzschild radius at the 
onset of this final 'delayed' collapse, the implosion will 
remain nearly spherical and generate little gravitational 
radiation. 

Now consider the alternative scenario, for which the 
thermal pressure generated by shock heating is insuffi- 
cient to support the star in quasi-static equilibrium. The 
star will then collapse promptly to a black hole on the 
relatively short dynamical timescale, and the bulk of the 
thermal energy and neutrinos will go down the hole. Un- 
able to achieve dynamical equilibrium, the compressed 
remnant will never relax to a spherical shape so that its 
implosion will be nonspherical. During the implosion, 
the mass may then emit a significant final pulse of grav- 
itational radiation before forming a static Schwarzschild 
black hole. 

We show below that the first scenario is energetically 
feasible, i.e, the merged remnant can reside in a hot, 
quasi-static stable equilibrium phase lasting many dy- 
namical timescales prior to cooling and contraction. Sur- 
prisingly, our conclusion is independent of the total rest 
mass of the object. Hence, even if the remnant mass 
greatly exceeds the maximum mass of a cold neutron star, 
the quasi-equilibrium will be stable and the collapse can 
be delayed considerably. 



correct equation of state of hot, nuclear matter is still un- 
certain, and the actual variation of T during compression 
and shock heating is not known. We note that T remains 
strictly constant and equal to 5/3 when ideal, degener- 
ate, nonrelativistic neutrons are heated to nonrclativis- 
tic temperatures by the shock. Following the passage of 
all the shock and acoustic waves and the accompanying 
generation of entropy throughout the matter, we approx- 
imate the function K(s) appearing in Eq. (Q) by a new 
constant, K(s) — if hot ||. We thus treat the merged 
remnant as a new, hot spherical polytrope once it set- 
tles down. Below we determine if hot > ifcoid uniquely in 
terms of -Kcoid for any n. 

Some of the above assumptions can be calibrated by 
comparing with existing Newtonian hydrodynamic simu- 
lations Q| [3|. The results of Rasio and Shapiro 
(hereafter RS; [2]), who treated the adiabatic case 
r = 2. are typical. They found that at most 5% of 
the total rest mass eventually escapes and that the rem- 
nant settles down to a hot, nearly spherical equilibrium 
state. The loss of energy via gravitational waves is about 
AE/M = 2AxlO- 3 M 1A R^ /2 , where M 1A is the mass of 
one star in units of 1.4M Q and i?io is the radius in units 
of 10 km. Hence, the assumptions of conserved rest-mass 
and total mass-energy during the initial dynamical epoch 
appear to be reasonable. 



III. KEY ASSUMPTIONS 



IV. NEWTONIAN ANALYSIS 



We assume that all of the original bulk kinetic energy 
at impact is dissipated into internal energy. Our goal is 
to assess whether the resulting pressure forces are capa- 
ble of holding a nonrotating spherical remnant in stable 
equilibrium. We treat the collision and merger adiabat- 
ically and we assume that the combined rest mass and 
total mass-energy of the system is strictly conserved. We 
neglect any net outflow of baryons from the system, as 
well as the loss of energy via gravitational waves and 
neutrinos. 

We adopt an adiabatic equation of state in the form 



P =K{s)p T Q , 



(1) 



where P is the pressure, po is the rest-mass density, 
r is the adiabatic index, assumed constant, and K is 
a function of the specific entropy, s. We assume that 
the approaching stars are cold with s — and take 
K(0) = i^coid to be a constant. We take the cold equation 
of state which supports the original neutron stars to be a 
polytrope, with polytropic index n given by T = 1 + l/n. 
We assume that the value of T defining the initial poly- 
tropic equilibrium and the value associated with subse- 
quent adiabatic changes are identical, as is assumed in 
most hydrodynamic numerical simulations of adiabatic 
gases. Adopting this assumption permits us to make a 
direct comparison with these simulations, as well as pro- 
vide a simplified analytic treatment of the collision. The 



When two identical polytropes, each of (rest) mass 
-^coid, radius i? C oid and gas constant -Kcoid, collide head- 
on from rest at infinity, the total initial energy of the 
system (excluding rest-mass energy) is given by 
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Here and throughout we set G = c = 1. The final rem- 
nant settles into a new polytropic configuration of the 
same index whose energy is given by 
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Conservation of rest-mass and total energy imply Mhot = 
2M co id and E^ot = -Ecoid , which, together with Eqs. (||) 
and (j^), require i?hot = 2i? co id- The mass-radius rela- 
tionship for a polytrope, 
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We test the validity of Eq. (§) by comparing with the 
Newtonian simulations of RS, who monitor the time- 
evolution of the total entropy for a T = 2 head-on col- 
lision (see Fig 7 of their paper). Once the remnant 
settles down, RS find that the entropy levels off at a 
value Sy(AicoidWM) = 21n(ii hot /A' cold ) sa 3. Hence 
Ahot/^coid ~ 4.5, which is in reasonable agreement with 
the analytic prediction of Eq. (jsj) , Ahot /^ cold « 4. 

V. GENERAL RELATIVISTIC ANALYSIS 

Here we reconsider the same problem in the context of 
general relativity. We again assume that the two poly- 
tropic stars collide head-on from rest at infinity, con- 
serving both their total rest mass, 2 Mo and total mass- 
energy, 2M. At t = the two spherical stars are widely 
separated so that their mass-energies, like their rest- 
masses, add linearly. We show that there exists a stable 
equilibrium solution for the remnant and determine the 
value of -Khot corresponding to this solution. 

Consider a sequence, parametrized by the central den- 
sity, of spherical equilibrium stars constructed from the 
same one-parameter equation of state. The equilibrium 
models are solutions of the Oppenheimer-Volkoff (OV) 
equation of hydrostatic equilibrium. Those solutions 
which are stable against radial collapse reside on the as- 
cending branch of the M vs. p c equilibrium curve, where 
dM/dp c > 0. Each of our colliding neutron stars is situ- 
ated on the stable branch of the OV spherical equilibrium 
curve for a polytrope of constant K — K co m at the start 
of its infall. We want to determine where the remnant 
resides on the corresponding curve for a polytrope of con- 
stant K = Khot once the stars merge and settle down (see 
Figl). 

The solutions exhibit a simple scaling behavior with 
respect to K: 

M(p c ) = K%m*(p* c ) , M (p c ) = K%m*(p* c ) 
R( Pc ) = K%r*(p* c ) , Pc = K- n p* c . (6) 

In Eq. (^) m*,mj and r* are universal, nondimensional 
functions of the nondimensional central density p* and 
depend only on the index n. The universal equilibrium 
curves for ann = l polytrope are shown in Fig 2. 

The rest-mass and total mass-energy curves have si- 
multaneous turning points at p* = /0* lax , at which point 
both masses take on their maximum values along the se- 
quence. The ratio m* jm^ measures the specific binding 
energy of a star according to -Ehind/Mn. = 1 — m*/m,Q. 
This universal ratio monotonically decreases with p* 
along the stable branch of the OV curve, p* < /5^ ax , 
as shown in Fig 3. Conservation of rest mass and total 
mass-energy imply 

M ,hot = 2M and M hot = 2M. (7) 

Combining Eqs. (||) and (0) gives 
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FIG. 1. Total mass-energy versus central density along the 
OV equilibrium curve for n — 1 polytropes. The solid line 
shows models for K = 1, the dotted line for K = 1.5 and 
the dashed line for K — 2. Solid points locate representative 
stable models with the same universal profile and binding 
energy. 
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FIG. 2. Mass versus central density along the universal OV 
equilibrium curve for an n = 1 polytrope. The solid line gives 
the total mass-energy, the dotted curve gives the rest mass. 
Solid points locate the representative stable model to which 
all the points in Fig 1 correspond. 
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FIG. 3. Binding energy parameter versus central density 
along the universal OV equilibrium curve for an n = 1 poly- 
trope. The solid point locates the representative stable model 
to which all three points in Fig 1 correspond. 



m 0,hot/ TO = m hot/ m * = 2 (#hcrt/-Kcold) 2 , (8) 

which implies 

KotMW = m*/m* . (9) 

The nondimensional binding energy ratio of the merged, 
hot remnant has the same value as in the cold progeni- 
tors. Because this ratio is a monotonic function of p* , the 
remnant can relax to the same stable, equilibrium point 
on the universal binding energy curve, hence 



l 0,hot 



ttiq and in. 



hot 



(10) 



Thus, up to a scale factor K, the remnant profile is iden- 
tical to that of its progenitors. Substituting Eq. ( |l0| ) back 
into m) implies 



^hot 



= 4" 



(11) 



whereby we recover the result obtained in the Newtonian 
analysis (cf. Eq. (||). 



VI. DISCUSSION AND CAVEATS 

Head-on collision and merger from rest at infinity need 
not lead to immediate collapse, regardless of the mass of 
the incident stars. Sufficient thermal pressure will always 
be generated to create a stable equilibrium configuration 
of higher entropy to which the fluid can can relax follow- 
ing merger: at fixed p 
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where P C oid is given by Eq. (pi) with K — K co \^ and where 



-fcold/Pcold ~ M/R. 



(13) 



Setting the thermal pressure in the hot remnant equal to 
the pressure of a Maxwell-Boltzmann neutron gas yields 
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(14) 



It is this thermal pressure that supports the remnant 
against collapse on a dynamical timescale (~ msec). The 
dominant cooling mechanism is thermal neutrino radia- 
tion, but the timescale is comparatively long (~ 10 sec) 
so the remnant evolves quasi-statically. 

It is straightforward to modify our argument to ac- 
count for the escape of baryons, the loss of energy due to 
gravitational radiation, the release of the progenitor stars 
from rest at finite separation rather than from infinity, 
etc. In the later two circumstances we simply replace 
Eq. (0) by 
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(15) 



where e < 1 measures the decrease in the total mass- 
energy below the initial value at infinite separation. Re- 
peating the argument provides a constraint on e for 
the remnant to relax to a dynamically stable quasi- 
equilibrium state: 
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where (m*/TOQ) m i n is the minimum value of the binding 
energy parameter along the stable branch. This mini- 
mum occurs at the maximum mass configuration and is 
« 0.9 for typical neutron stars (see Fig. 2). The low 
expected radiation efficiency / = 1 — e<l jjj j| || 
implies that all head-on collisions from rest at infinity 
can lead to quasi-equilibrium except for those progeni- 
tors with masses within / of the maximum. For colli- 
sions starting from rest at finite separation d, the total 
mass-energy is decreased by the initial binary interaction 
energy, so that 



hot. 



2M -M 2 /d, 



(d > M), 



(17) 



Equating Eqs. ( [l5|) and ( |r7| ) and ignoring radiation 
losses implies e w 1 — M/(2d), whereby the condition 
for stable equilibrium becomes 



d > 



M 



m 



(18) 



I 1.2) As M M max and (m* /ml) — > (m*^TOo) min , we require 



that d — > oo for stable equilibrium 
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Simulations of head-on collisions of neutron stars in 
full general relativity have been performed || @, but 
no systematic studies have been published. The reported 
cases exhibit dynamical behavior qualitatively similar to 
the Newtonian simulations regarding the recoil shocks, 
bounce and damped oscillations following the collision, 
and there is no collapse This result is consistent 

with our theoretical expectations, but we await a more 
complete numerical investigation for confirmation. There 
is the possibility that the configuration might greatly 
overshoot the final, allowed equilibrium state during the 
merger process and plunge immediately to a black hole, 
but such a fate is not evident. Hopefully, this paper may 
stimulate further numerical experiments, including ex- 
amples which employ a more realistic hot, nuclear equa- 
tion of state. Off-axis collisions are also of interest; while 
shock heating will be less important in this case, the an- 
gular momentum acquired by the remnant may supple- 
ment the reduced thermal pressure to prevent sudden col- 
lapse |^] . Resolving the fate of a head-on collision might 
serve as a proving ground for a relativistic code designed 
to handle these more complicated 3D binary scenarios. 



[8] The gravitational waveform for the relativistic simulation 
reported in Refs ]B| matches well with the Newtonian 
result reported in Ref g ; see Fig 8 of Ref Q . The wave- 
form probes the bulk motion of the coalescing matter and 
therefore attests to the similarity of the behavior. 

[9] T.W. B aumgarte and S.L . Shapiro, Astrophys. J. in press 
(1998); |astro-ph/9801294| 
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